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Abstract 

The High Density Effective Theory formalism is used to calculate 
the low energy properties of the phonons and gluons in the Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) phase of two-flavor QCD. The ef- 
fective phonon Lagrangian for the cubic crystal structure, which is 
favored according to a recent study, depends, at the second order in 
the derivatives, on three parameters which we calculate in this paper. 
We also compute for soft momenta the effective lagrangian for the 
gluons of the unbroken SU(2) C group, both for a single plane wave 
structure and for the cube. In both cases the Meissner gluon mass 
vanishes as in the case of complete isotropy; on the other hand there 
is a partial Debye screening due to the existence of blocking regions 
on the Fermi spheres. For the single plane wave structure the gluon 
dielectric tensor is non isotropic, while it is isotropic for the cubic 
crystal, in spite of the intrinsic anisotropy of the structure. 

1 Introduction 

The study of the color superconducting phases of QCD |L] has recently 
received much attention and is at the moment one of the most studied frontier 
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problem of QCD ||. Cold QCD with three degenerate flavors, for very large 
values of the quark number chemical potential /x, is in the so-called color- 
flavor locked (CFL) phase J4|. In this phase, pairs of quarks of opposite 
momenta, each of them lying near the Fermi surface, condense according to 
the BCS pairing scheme known from ordinary superconductivity. 

The CFL phase is expected to persist also below asymptotic values of 
/i, when the flavor degeneracy is broken, with the strange quark having a 
larger mass than the lighter quarks. In absence of electrons, the CFL phase 
is electrically neutral. It has the same number of u, d, and s quarks, which 
corresponds to the lowest energy for the Cooper pairs. 

When the value of \i is further decreased, one expects that pairing between 
quarks of non vanishing total momentum becomes possible. Each quark 
momentum will lie near its Fermi surface (Fermi spheres will now be of 
different radii) [|5], |5], |7|, K |S| . Closeness to their Fermi surfaces for both quarks 
in the pair will correspond to a lower free energy. Such a possibility for the 
quark pairing leads to the suggestion of a crystalline color superconductivity 
H H 0' B HI- We assume the existence of this phase and we will limit our 
analysis, as in the existing literature, to two flavors, with 5fi 7^ (5{i the 
difference in the quark chemical potentials). 

In certain condensed matter systems (ferromagnetic systems with para- 
magnetic impurities) [TI| the possibility of fermion pairing of non vanishing 
total momentum had already been discussed. The name LOFF phase is 
therefore commonly used, formed from the initials of the authors of ref. [TO 



Translational and rotational invariance are spontaneously broken within 
a LOFF phase. The gaps are expected to vary according to some crystalline 
structure. The determination of such a structure is obviously of the upper- 
most importance and this problem is discussed and approximately solved 
by Bowers and Rajagopal [O] in a recent paper. To understand the crys- 
talline superconductivity one has indeed to know the space dependence of 
the Cooper pair condensate. A simple plane wave ansatz had been used for 
a preliminary understanding of the possible dynamics. The plane wave has 
the finite momentum of the Cooper pairs, on the assumption that all the 
possible Cooper pair momenta choose a unique direction. These studies (see 
H and (TOJ) already lead to the conclusion that LOFF pairing is favored, 
over BCS pairing or absence of pairing, within a certain window in 5fi. The 
window starts with a first order phase transition to the crystalline phase. It 
ends with a second order transition where translation invariance is restored. 
The window for QCD with two flavors had been found in ref. || by modeling 
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quark interactions via four-linear Fermi couplings among quarks and for the 
case of a plane-wave condensate structure. A much wider window is found 
however in a model based on gluon exchange ||. This follows essentially 
from the fact that gluon exchange favors forward scattering, rather than the 
s-wave scattering typical of the Fermi couplings. The effective low energy 
phonon lagrangian for the case of the single plane wave structure has been 
discussed in [J/J. Since translation invariance is only partially broken one 
phonon is present. 

A comprehensive study of the different possible crystal structures for the 
crystalline color superconducting phase of QCD is presented in JTT]] where 
the authors perform a Ginzburg-Landau expansion in the vicinity of the 
phase transition. Instead of the single plane wave ansatz they now introduce 
a multiple plane-wave ansatz and find that a face-centered cubic crystal is 
favored. This structure is built up by eight plane waves directed to the 
corners of a cube. At least qualitatively it seems true that the cubic crystal 
remains the favored structure also in complete QCD (although not perhaps at 
very asymptotic densities). It is important to note that for the cubic crystal 
structure three phonons should be present since translational invariance is 
completely broken. 

The crystalline color superconductivity has a relevant place in the QCD 
phase diagram and as such it deserves to be studied and to be known in 
its main features. Moreover it can play a role in the physics of compact 
stellar systems where, within the core, a crystalline phase might be reached. 
The authors of ref. |J had already pointed out that, if the compact star 
is rotating, the ensuing rotational vortices may become pinned and as such 
play a role in the formation of glitches. This subject is certainly interesting, 
but we shall not discuss it here (for a review see ref. ]12] ) . 

We have already mentioned that crystalline superconductivity had ini- 
tially been proposed in condensed matter physics. The QCD studies, in 
particular also the present one, may be translated into corresponding prop- 
erties for condensed systems, at least for the phonon properties. We mention 
in particular the field of atomic systems, gases of fermionic atoms below their 
Fermi temperature ]13], An advantage of testing the ideas of crystalline 



superconductivity on these systems is that their parameter space is less rigid 



than for QCD. For example see ref. |15 | 



The plan of the present work is as follows. In section 2 we first briefly 
review the formalism we will use in this paper, i.e. the High Density Effective 
Theory (HDET). This effective theory is based on the limit —>■ oo and uses 
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velocity dependent fields in the same spirit of the Heavy Quark Effective 
Theory (HQET) |L6| where the limit rriq — > oo is taken. It is discussed 
in [H], [18], |19| and more extensively in ref. f20fl , where the method and 



the notations are explained in more detail. In section 3 we discuss the low 
energy dynamics of the phonons. In section 4 we calculate the parameters of 
the effective phonon lagrangian. In section 5 we study the properties of the 
gluons of the unbroken SU(2) C for the case of the plane wave structure and 
in section 6 we perform the same calculation for the LOFF condensate with 
the cubic symmetry. The application to the cubic crystal structure requires 
some group theoretical definitions which are summarized in Appendix 1. In 
Appendix 2 and 3 we report some details of the calculations. 

2 Heavy Density Effective Theory in the LOFF 
phase 

We begin with the effective lagrangian of HDET for gluons and ungapped 
massless quarks. If Cq is the free quark lagrangian and L\ represents the 
coupling of quarks to one gluon, the high density effective lagrangian can be 
written as 

C D = C + C 1 + C 2 + {L^R). (1) 
Let us begin with the first two terms, i.e. 

C = ty+iV ■ dtp + + ipliV ■ dtp.] , (2) 
C x = i g J2 [p\iV ■ Aip + + ipliV ■ A^j_] . (3) 

v 

The velocity dependent left-handed field ip + = ipv is the positive energy part 
in the decomposition 

Hx) = J ^e~^ x lMx) + M*)\ ( 4 ) 
while tytf is the negative energy part which has been integrated out. Here 

= e^P+^x) = J ^ e -^P + m (5) 
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and 

9«(x) = e^P.^x) = J ^e-^P.m . (6) 



P± are projectors defined by 



P± = \{l±a-v) . (7) 

The cut-off 5 satisfies 6 <C fi while being much larger that the energy gap. 
These equations correspond to the following decomposition of the quark mo- 
mentum: 

p K = ^ + t (8) 

where 

= (0, v) (9) 

and \v\ = 1. We also use V 1 = (1, v) ,V" M = (1, —v). 

The other velocity-dependent field ip- = ip-v is obtained from ip + by the 
change v — ► — {£. Therefore the average over the Fermi velocities is defined 
as follows: 

?-/£■ (1 °> 

The extra factor - is introduced because, after the introduction of the field 
2 

with opposite velocity ■?/>_ , one doubles the degrees of freedom, which implies 
that the integration is only over half solid angle. 

The effective lagrangian (|J) is obtained in the limit \x — > oo [17, 18, 19, 



PU[ . Therefore Cq and C\ do not contain the chemical potential because the 
dependence on \x has been extracted from the quark field operators. We also 
note that one has a single sum over the velocities, because the quark fields 
are evaluated at the same Fermi velocity; as a matter of fact the off-diagonal 
terms are cancelled by the rapid oscillations of the exponential factor in the 
\i — > oo limit. This point will be discussed in a more detailed way in the 
sequel. 

Besides Cq and C\ eq.([T]) contains the term 



C 2 = -Y P^ v W+ = D u D u ib + + ibl D u D u ib_ 



C 2 is a non local lagrangian arising when one integrates over the ^ degrees 
of freedom in the functional integration. It contains couplings of two quarks 
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to any number of gluons and gives contribution to the gluon Meissner mass 
(see below). We have put 

pf iu _ gl xv _ + yvytA^ _ ( 12 ) 

This construction is valid for any theory describing massless fermions 
at high density provided one excludes degrees of freedom far from the Fermi 
surface. Let us now examine the modification in the formalism in the presence 
of a LOFF condensate; to begin with we consider the approximation of one 
plane wave. The space dependence of the crystalline condensate in this case 
is 

A(f) = A exp{2ig • x} . (13) 

The effect of the non vanishing vacuum expectation value can be taken into 
account by adding to the lagrangian the term: 

£ A = -— exp{2iq ■ x} e a p Z eijil^ a {x)Cij)jp{x) - (L -»• R) + h.c. . (14) 

In order to introduce velocity dependent positive energy fields Vfyjia with 
flavor i we decompose both fermion momenta as in (§) and we have: 

£a = — g ex P{^ ■ a(vi, vj, ^}e ij e ot ^'4F^. ia {x)C'4)-^. j p{x) 

Vi,Vj 

—{L —> R) + h.c. , (15) 

where 

a(vi, vj, q) = 2q- - fijVj . (16) 

Let us define 

H= 2 ' ^ = 2 ' ^ ' 

We will consider large values of /i (with 5fi/ fi — > 0, q/Sfj, — > and g fixed), 
but below the value where the CFL phase sets in. To give our procedure a 
mathematical meaning, we consider a smeared amplitude as follows: 

lim expizr • a(vu v 2 , q)\ = lim / df" expizr" • a(vu v 2 , q)}q(r,f') . (18) 

fl^OO fl— >oo J 

In [^TJ we assumed the following smearing function: 

3 

g(f,f') = gi(r- r) = JJ 6 M {r k - r' k ) (19) 

k=l 



where 



e M (r - r') 



1 



9 r' — r + 



At 



9 r' 



At 



r — 



_i_ f or | r ' _ r ] < A€ 



(20) 







elsewhere 



(6* the Heaviside function). In the present paper, for reasons to be discussed 
below, we will choose 



sin 



-Kq{r k ~ r' k ) 
R 



(21) 



We evaluate (|18|) in the /i — > oo limit by taking <f along the z— axis, and 
using the following identity: 



J d 3 r" exp{if" ■ a}gjj(r — f') = exp{if ■ a] ^— J <5fj. 



(22) 



where 



8r{x) 



R r I T 

— tor x < — - 

7T ' ' 2i? 



(23) 



elsewhere 
For the components x and y of a we get 

nq 



KfJLlVi + H2V2)x,y\ < 



R 



i.e. approximately (for 5/x <C ji) 



\{Vl + V 2 ) x ,y\ < — . 

From this in the high density limit it follows 

iTi = -iT a + 0(1///) . 



(24) 



(25) 



(26) 
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More precisely, if 9\ and 62 are the angles of v% and with respect to the 
z-axis one gets 

1 = 2 +7T + 2— tan0 2 - (27) 
For the z component we get 

a z = 2qh{cos6 2 ) . (28) 



The expression for h(x) is [21]: 



+ (29) 



Neglecting corrections of order Sfi/fj, one finds 

h < x ) = 1-^*1, cos e= 6 -^. (30) 

x q 

The two factors tt/R arising from the x and y components are absorbed 
into a wave function renormalization of the quark fields, both in the kinetic 
and in the gap terms. As for the z component one remains with the factor 

^ hz 5 R [h(v-n)) » ^S R [h(v- ft)} (31) 

where 

n = q/q (32) 

in the gap term, whereas for the kinetic term we get a factor of 1. 

As we noted already, in ref. [[21] we used as a smearing function gj (see 
eq. (|i9|)) and 5r turned out to be given by 

5 A (x) = ^^, R = q\A£\. (33) 

Both eqns. (|33|) and (|23|) define regions where 8r 7^ 0, i.e. define domains 
where pairing between the two quarks can occur; they correspond to the 
pairing region (in contrast with the blocking region, where Sr = 0), using the 
terminology of |4j]. However, given the results of [|ll|], which show that the 
different sections of the pairing domain for the LOFF cubic structure are non- 
overlapping, we prefer to choose the prescription (|2l|) and, consequently, eq. 
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for Sr. We explicitly observe that, owing to (PUD, with x = v ■ n = cos 9, 
eq. fl23|) defines a 'ring' on the Fermi sphere. Our approximation consists in 
putting exp[i2qhz] = 1 in eq. (|3T|), owing to the presence of 5r that enhances 
the domain of integration where h = 0, for R sufficiently large. 

Let us discuss the range of values of i?Q for the plane wave ansatz in more 
detail. We notice that R = oo implies the vanishing of the pairing region. 
The precise value of R could be fixed by the gap equation. For instance, 
by using the fact that near the second order phase transition the pairing 
region has an angular width of order A/5fi (see for a discussion |Tl|), from 
the definition of 5r we find (for A — > 0) 



7T _ A 

R 5fi\\5[i 




1 . (34) 



This shows that R — > oo for A — > 0. In any case for the purpose of this 
paper we leave R clS cL parameter. 
In conclusion 

A 7T 

£ A = ~— ^5 R [h(v-n)}e ij e a ^i>l ia (x)C'ilj_^ jl3 (x)-(L -> i?) + h.c. (35) 
We now change the basis for the fermion fields by writing 
V> +1 « = E^P^ (*, « = 1,2) 

^+,31 = 

^+,32 = 4>l, (36) 

where are the Pauli matrices for A = 1, 2, 3 and Co = 1; similar expressions 
hold for ipj:. 

A different, but also useful notation for the fields 4>+,ai uses a combination 
of A matrices, as follows 



^ = E • (37) 

A=0 V 1 

The matrices are defined in terms of the usual A matrices as follows: 
A = ^A 8 + y|Ao; Xa = \a(A= 1, 2, 3); A 4 = A 5 = In this 

1 Notice that the quantity R used in the present work is not the same used in pjj. 
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basis one gets 



Here 



= E E x 

v A,B=0 



At ( iTr[f\ V ■ D f B ] 
A AB 



A 



AB 



iTr[T\V ■ D* T B ] 



X 



V2 VC^. 



and 



r A = ^ (A = o,...,5). 



J* 

(38) 
(39) 

(40) 



where the matrix A^b is as follows: A AB — (AorS = 4 or 5), and, for 
A,B = 0,...,3: 



A 



AB 



A 



(41) 



having put 



( 1 ^ 
0-10 
0-10 

V o o o -l / 

A eff = ^5 R [h{v-n)]. (42) 



AB 



From these equations one can derive the propagator for gapped quarks: 

1 



D AB {£) 



V-£V-£ 



A 2 



V-£6 AB -A AB 
-A AB V-£S AB 



(43) 



On the other hand the propagator for the fields x 4 ' 5 does not contain gap 
mass terms and is given by 



D(£) 



(v-iy 1 
o (v-e)- 1 



(44) 



These formulae can be easily generalized to the case of the face centered 
cube |Ll||. This is characterized by eight qj vectors with equal length q = n/a 
and directions corresponding to the eight vertices of a cube. 

In Fig. 1 some of the symmetry axes of this cube are shown: they are 
denoted as C 4 (the three 4-fold axes), C 3 (the four 3- fold axes) and C 2 (the 
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six 2- fold axes). For more details on the symmetry group of the cube see 
Appendix 1 and ||22|| . The space dependence of the condensate corresponding 
to this lattice is as follows 

8 _ ( 2iri 1 
A(x) = A Jj exp{2iq k -x} = A ^ exp < (dXi + e 2 x 2 + e 3 x 3 ) \ , (45) 

fc=l e i= ± L a > 

where the eight qj vectors with equal length 

q = n/a (46) 

have directions corresponding to the eight vertices of a cube. 




Fig. 1 - The figure shows some of the symmetry axes of the cube, denoted as 
C 2 , C 3 and C 4 (see the appendix for details). We show also the parameteri- 
zation used for the coordinates of the vertices. 



To describe the quark condensate we add the term: 



A 8 

C A = -— exp{2iq k ■ x} e a pzeij^J a {x)Cil)jp{x) - (L 
1 k=i 



R) + h.c. , (47) 



completely analogous to (|T4T). 
Let us put 



nk = qk/\qk\- 



(48) 
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By the same procedure used for the plane wave condensate one has 



A 8 7T 

£a = -— ^ S R[K^-^k)]eije a ^^. ia (x)Cil)-e ; j0{x) 

-(L^i?)+h.c. (49) 

In ( |48D the eight unit vectors defining the vertices of the cube are 
n 3 = -^=(-l,-l,+l), n 4 = -L(-i,+i,+i), 

TJr = i(-l,-l,-l), ffe = i(-l,+l,-l). (50) 
C + L\ + £a is still given by eqns. (p8|)-(p|) but now 

^e ff = ^jZHKv-n k )} ; (51) 
K k=i 

the quark propagator is given by Q43Q with A e jj given by (|5T|) . 

An interesting point can be noted in connection with eq. fl51"D. This 
equation shows that the pairing region for the cubic LOFF condensate is 
formed by eight distinct rings; each ring is associated to one vertex of the 
cube and has as its symmetry axis one of the C3 axes of fig. 1. According to 
the analysis of the LOFF vacuum state corresponds to a situation where 
these domains have at most one common point. Given the symmetry of the 
cubic structure we can limit the analysis to one pair of rings, for example 
those associated to the vertices Hi, in fig. 1. The common point between 

these two rings lies on the axis C 2 of fig. 1 and has v = — ^=(1, 1,0). Since 

v2 

it must also belong to the boundary of the two pairing regions we have the 
condition (p3|): 

\h(v-n 1 )\ = ^, (52) 

which implies 

R = . (53) 

2h[J2/3 
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In this particular case the approximate form for h(x) (see eq. (|3T)|)) is not 
appropriate owing a numerical cancellation. In fact for q/\8[i\ ~ 1.2, \x = 
400 MeV, 5li = 40 MeV, one would get R « 76, whereas the use of the 
complete equation (p9| ) gives 

i? ss 18 . (54) 



3 Phonons in a cubic LOFF crystal 



The condensate fl4"5|) breaks both translations and rotations. It is however 
invariant under the discrete group Oh, the symmetry group of the cube. This 
can be seen by noticing that the condensate is invariant under the following 
coordinate transformations 



R x : 




R 2 : 


Xi 


R 3 : 


Xi 


I : 





-x 3 , X 2 

%2, X 2 — 
— Xi, x 2 



x 3 , x 3 -> -x 2 , 
-> x 2 , x 3 —± X\, 
X\j x 3 > x 3 , 
-> -x 2 , x 3 -> - 



(55) 



that is rotations of 7r/2 around the coordinate axes, and the inversion with 
respect to the origin. Since the group Oh is generated by the previous 4 
elements (see Appendix 1) the invariance follows at once. 

The crystal defined by the condensate ( fig) can fluctuate and its local 
deformations define three phonon fields <f>® that are the Nambu-Goldstone 
bosons associated to the breaking of the translational symmetry. They can 



be formally introduced by the substitution in (|45|) 

2n 



2qx i -> $ (i) (x) = — 
a 

where the three scalar fields $W satisfy 

($ W )o = 

whereas for the phonon fields one has 

<0 W (*)>O 



x l + 



*)//• 



2vr 



- ,1 7 ; 



. 



(56) 



(57) 



(58) 
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In ( |5*6] ) / is the decay constant of the phonon with dimension of an energy. 

di) 

k\ k2 K3 



We have therefore three fluctuating fields <pi, k „u„ for any elementary cube 



defined by discrete coordinates 

k\TX k 2 7T fc 3 7T 

a^fci = , Vk 2 = , Zk 3 = , (59) 

q - q q 

i.e. 

4>k ) 1 k 2 k 3 = (p (l) (t,x kl ,y k2 ,z k3 ) . (60) 
The interaction term with the NGB fields will be therefore given by 

r /ttX 3 +°° +°° +°° 8 , . 

Sint = -Jdtl-) E E E EEMi^'j/} 

V " / ki=— oo k2=— oo kz=— oo v m=l 

x eye^VLff^^-if- + /i.e. (61) 

where 

3 

k2 kz Z-^i i rki k2 kz 

(62) 

»=i 

and the eight vectors e^" 1 -* are given by 

(ef } ) = e M = ^- (63) 

Notice that translating the fields in momentum space one would have to con- 
sider Fourier series with respect to the discrete coordinates. This introduces 
in the theory a quasi-momentum p through the combination exp(i(p x x kl + 
PyVk 2 + VzZkz)- Therefore one has to restrict p to the first Brillouin zone 
-q < Px < q, ~q <p y < q, -q < Pz < q- 

The complete effective action for the NGB fields <f>® will be of the form 



. / \ 3 -(-oo -(-oo -(-oo 

S=ldt[-\ E E E £(<t> {i) (t, hir/q, hTi/q, k^/q)) . (64) 

V ^ / fci=— oo k2=— oo &3=— oo 



In the action bilinear terms of the type 4> klk2k3 4>k\ k ' 2 k' 3 w ith {fci &2 ^3} 7^ 
fcj fc 2 fc 3 may arise. In the continuum limit this terms would correspond to 
partial derivatives with respect to the three space directions. 

All these considerations become superfluous when one realizes that we are 
actually interested in an effective description of the fields (j) kl k2 ks in the low 
energy limit, i.e. for wavelengths much longer than the lattice spacing ~ 1/q. 
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In this limit the fields 4> ki k2 ks vary almost continuously and can be imagined 
as continuous functions of three space variables x, y and z. Therefore we 
shall use in the sequel the continuous notation r) . 

The coupling of the quark fields to the NGB fields generated by the 
condensate will be written as 

A^ T C^ ex P {*(ei$ (1) + e2$ (2) + e 3 $ (3) )} (65) 

e<=± 

making the theory invariant under translations and rotations. These invari- 
ances are broken spontaneously in the vacuum defined by eq. (|57|). 

Integrating out the Fermi degrees of freedom one gets an effective la- 
grangian containing only the light degrees of freedom, i.e. the gluons of 
the unbroken subgroup SU(2) C and the phonons. We will discuss gluons in 
Section ^| and we consider here only the phonon effective lagrangian. 

The effective lagrangian for the fields $W has to enjoy the following sym- 
metries: rotational and translational invariance; Oh symmetry on the fields 
$W. The latter requirement follows from the invariance of the coupling ( |65| ) 
under the group Oh acting upon $w. The phonon fields 4>^\x) and the 
coordinates x % must transform under the diagonal discrete group obtained 
from the direct product of the rotation group acting over the coordinates 
and the Oh group acting over $W(a;). This is indeed the symmetry left after 
the breaking of translational and rotational invariance. 

In order to build up the effective lagrangian, we start noticing that if X^ 
are quantities transforming as the same representation of the $W's under Oh 
(a 3-dimensional representation), then all the invariant expressions can be 
obtained by the following three basic invariant expressions 

— ^i+-^2+^3' ^4pQ) — X X X 2 +X%X$ +X 3 Xl, Ie(Xi) = X^X^X^, 

(66) 

which is a general property of the symmetric functions of three variables. The 
invariance is easily checked by noticing that these expressions are invariant 
under the three elementary rotations Ri and the inversion / (see eq. (|55|) ). 
Quantities transforming as 3>W are 

(JjL 

We proceed now as in ref. 0, by noticing that the most general low-energy 
effective lagrangian, invariant under rotations, translations and Oh is given 
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a) b) 
Fig. 2. Self-energy (a) and tadpole (b) diagrams. 



by 

£ = C E (^ W ) 2 + A(/2(V$«), / 4 (V$«), J 6 (V$«)) . (68) 

^ 1=1,2,3 

The reason for allowing all the possible spatial gradients in this expansion is 
that the vacuum expectation value of V$^- ) is proportional to q and it is not 
in general small. By following ref. |23| one gets the effective lagrangian for 
the fields 0^ at the second order in the derivatives as follows 

c = \ E (0 W ) 2 -^ E IV0«| 2 -^ E (d^r-c e d^%^. 

i=l,2,3 i=l,2,3 i=l,2,3 i<j'=l,2,3 

(69) 

The effective theory depends on three arbitrary parameters, a result that 
follows also from the general theory of the propagation of elastic waves in 
crystals P3 . 



4 Parameters of the effective lagrangian for 
the phonon fields 

In this Section we will compute the parameters a, b, c appearing in (p5|). We 
use the gradient expansion. First we introduce external fields with the same 
quantum numbers of the NGB's and then we perform a derivative expansion 
of the generating functional. This gives rise to the effective action for the 
NGB. At the lowest order one has to consider the diagrams in Fig. 2, i.e. 
the self-energy, Fig. 2a, and the tadpole, Fig. 2b. The diagrams can be 
computed observing that, by expanding (|65D in the phonon fields, one has 
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couplings of two quarks with one and two phonons. Let us put 

3 



^(t,T0 = E^V (i) (*.r)- 



(70) 



i=l 



At the first order in the fields one gets the coupling with one NGB 



ix R 

(L^R) + h.c. 



f 



(71) 



and the coupling of two quarks and two NGB's 
C 



8 7rA (cp( m )) 2 
H m=l U Z J 



- (L-yR) + /i.e. 
In the basis of the x fields: 



gt + g\ 



(72) 



(73) 



Here 



#3 



#4 



8 7rA 
E -5-^fl[M^-"m)]— 
m=l K J 



E -5" <W( U • n m)\ of2 
m=l K 2 J 



{ 1 ^ 

0-10 
0-10 

V -1 ) 

{ 1 \ 
0-10 

0-10 

V -1 J 



AB 



(74) 



AB 



To perform the calculation one employs the propagator given in Eq. (p!3|) 
with A e jj given in (|51| ) and the interaction vertices in (0). The result of 
the calculation of the two diagrams at the second order in the momentum 
expansion is as follows: 



* v m,k=l 
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x 8 R [h(v-n k )](itp {ky 



x 



D(£)D(£ + p) 
2A 2 eff + V ■ £V ■ (£ + p) + V ■ £V ■ (£ + p) 



r ( \ ■ 4 x Vv- ^ /" d 2 l 7rAA e// _ (m) 2 

£eff(p)tad = ij^rpY, Ey ^ - •"*>)] j ) 

where 

D(£) = ^-4-A e 2 // + ^e , (75) 

One can easily control that the Goldstone theorem is satisfied and the phonons 
are massless. As a matter of fact one has 

. 4 x 4/i 2 vrA 

Cmass = C e ff{U)s.e. + £>eff{u)tad = l ^3^2 X 
r rl 2 f r 8 7rA 

t? V / m,k=l 

+ A e// E^[^(^-«m)](^ (m) ) 2 ] ■ (76) 

m=l 

Both terms in the r.h.s. of eq. ( |76[ ) present a double sum: m, = 1, 8 
(see eq. fl5lD for A e jf); we have seen at the end of section || that the eight 
rings that define the pairing region are non overlapping. Therefore only 
the terms with m = k survive and one immediately verifies the validity of 
the Goldstone's theorem, i.e. the vanishing of ([76]) . Notice that in our 
approximation the masses of the Goldstone bosons vanish only if the pairing 
regions are not overlapping, signaling that when the regions overlap we are 
not at the minimum of the free-energy (see [p"T|j). 

At the second order in the momentum expansion one has 

£e ff (p) = *^r^E E \ Cf) 2 HKv-n m m^) 

•> v m,k=l 

x ^.^(i^y'^ g^^f . ( 77 ) 

Using the result 

iT (78) 
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and the absence of off-diagonal terms in the double sum, we get the effective 
lagrangian in the form 



£eff(p) 



7T 



A x2 



(fc) 



x (v-p) ip {k) (y-p)<p 



(6 R [h(v ■ n k )]f 
A 2 



(79) 



To perform the calculation one can exploit the large value found per R , i.e. 
eq. fl54]) . Since for R — > oo the <5r function becomes the Dirac delta, we take 
the .R — > oo limit and we handle the function according to the Fermi trick 
in the Golden Rule; in the numerator we substitute one 5 R function with the 
Dirac delta while for the other one we take 



n6 R [h{x)] tt5 r (0) 



R 



R 



(80) 



We observe that this limit is justified by the large value found for R, see eq. 
fl54|). The sum over k in (ff9|) gives 

E 6 R [h(0-n h )]<pW6 R [h(0.fi k )]<pW -> § E s[h(v-n k )} (v {k) ) 2 = 



k=l 



71 



k=l 



R 



71 



E 6 

k=l 



5ji 



qv ■ n k 



(k) 



with 



k R 



2 R< 



7P 



fc=l 



v ■ n k - 



5p 



^))-(81) 



qR 



(82) 



On the other hand A^, in the denominator of (|79|) gives 



A 2 

eff 



A . 



(83) 



Therefore one gets 

..27, 3 8 

i,j=l k=l v 



£>eff(p) 



2vr 2 / 2 ^ 



q 



Vj/J^efp^p^ 



(84) 

The integration over the Fermi velocities require special attention. We use 
the result 



J2 e (k) e (k) 



~i 3 



86 



(85) 



k=l 
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this fixes the constant multiplying the time derivative term in the effective 
lagrangian at the value (taking into account (|10D) 



R 



4vr 2 / 2 



(86) 



Therefore one obtains canonical normalization for the kinetic term provided 

2 8n 2 k R 

f ~ ~2^~ ■ (8?) 

To compute the parameters a, b, c of the effective lagrangian (|69| ) we need 
to evaluate 



8 



^'^/^'P,} ' (88) 

We find the result (see Appendix 2) 

0L = jjWy + 3COS2 3 ^~ 1 £ ^ (89) 

where 

Pij = $is8jt + 8 it 5 js , r, s, t in cyclic order . (90) 
From this equation we get 

1 / n2 ,1*2 1 



£«//(?) = 2 (p°V (,)2 - g «yp m ^"v u> ~ 

;(^*-|>'- w ,f'~ 1 £ pVW)- < 91 > 

i<j=l,2,3 



Comparing with eq. (|69"D we finally obtain 

1 , n 3cos 2 #„-l 
°=12- " = - C= iJ— ■ (92) 

One can note that these values for the parameters of the phonon effective 
lagrangian satisfy the positivity constraints found in 



In p3| one can find the general discussion of the dispersion relation for 
arbitrary values of the parameters a, b and c. In the present paper we give the 
expression for the eigenvalues of E 2 , according to the direction of propagation 
of the momentum p, in table 1. 
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Momentum 


Eigenvalues 


Eigenvectors 


c A 


p = pe w 


1/12 
1/12 


0» 


c 3 


o = p(e( 1 ) + e( 2 ) + e (3) ) 


(1 + 3 sin 2 0,)/6 
(1 + 6 cos 2 9 )/36 


0(1) + 0(2) + 0(3) 




(1 + 3 sin 2 0,)/6 
(1 + 6 cos 2 0„)/36 


0(D +a; 0(2) + ( 1 _ ;r )0(3) 
_0(1) + 0(2) + 0(3) 


p = p(e«-e (2) + e (3) ) 


(1 + 3 sin 2 0,)/6 
(1 + 6 cos 2 9 )/36 


0(1) +X( f)(2) _ (l_ x )0(3) 
0(1) _ 0(2) + 0(3) 


p = p(e( 1 ) + e (2) -e (3) ) 


(1 + 3 sin 2 9 q )/6 
(1 + 6 cos 2 9 q )/36 


<f)W + X(j)W + (1 + x)0( 3 ) 

0(1) + 0(2) _ 0(3) 


c 2 


p = p(e (i) + e (? ' ) ) 
i < j 


1/12 
sin 2 9 q /8 
(1 + 3 cos 2 9 )/24 


0« 

+ 

0(0 + 0(i) 


p = p(e {%) -e U) ) 
% < j 


1/12 
sin 2 9 q /8 
(1 + 3 cos 2 ? )/24 


0« 
0« + 00) 
_0« + 0(i) 



Table 1 - Eigenvalues and eigenvectors of t-^; p along the symmetry axes. 



5 Gluon dynamics in the the LOFF phase: 
single plane wave structure 

In this Section and in the subsequent one we wish to derive the effective 
lagragian for the gluons of the unbroken SU (2) c subgroup of the two-flavor 
LOFF phase. To start with we review the results obtained for the gluon 
effective lagrangian at one- loop in the 2SC isotropic phase [||. They were 
obtained in p5| ; in |26j the results of J25] were obtained by the method 
of the High Density Effective Theory and extended to the three-gluon and 
four-gluon vertices. In the 2SC phase one considers massless quarks of two 
different flavors at the same chemical potential /i and at low temperature. 
Color superconductivity breaks the color group SU(3) C down to a SU{2) C 
subgroup and the quarks with nontrivial SU{2) C charges acquire a gap A . 
Therefore only the three gluons relative to the generators of the unbroken 
subgroup remain massless. 

The lagrangian including one loop corrections for the SU(2) C gluons can 
be written as |2§, |2§] 



C = ^(E?E? - flfflf) + ^E?E? , (93) 
with a = 1, 2, 3 color indices and 

k=^L- 04) 

Eqns. (|93|) and fl94|) mean that, at one loop, the medium has a dielectric 
constant 

e = k + l (95) 
and a magnetic permeability A = 1, while the gluon speed in the medium is 

' oc^. (96) 



As noted in [pSj , it is possible to put (^) in a gauge invariant form 

C = -\F^F; V (a = 1,2,3), (97) 

provided the following rescaling is used 

A a -> Aq = k 3 / 4 A a , (98) 
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A* -> Af = k 1/4 A1 , (99) 

xo x' Q = h~ l l 2 XQ , (100) 

g^g' = k-^g . (101) 

Let us now consider the LOFF non-isotropic phase. To begin with we 



assume the crystal structure given by a plane wave i.e. A(x) given by 
The difference with the isotropic case consists, in the present case, in the 
substitution 

A — »• A ef f = ^-S R [h(v ■ n)] , (102) 

where n = q/q and h{x) given by (|30| ). 

The effective action allows the evaluation of the one loop diagrams with 
two external gluon lines and internal quark lines similar to those in fig. 2. 
Their contribution to the polarization tensor is written down in Appendix 3. 

Let us now write 

iC(p) = iC(o) + *iC(p) . (los) 

As discussed in Appendix 3 we find a vanishing Meissner mass 

IL? 6 (0) = (104) 
and a non- vanishing Debye screening mass 



gfi cos d a - cos 9 b 

m D = — yi + 2 ' ( 105 ) 

where cos^ a and cos 9b (—1 < cos^ a < cos^ < 1) are the solutions of the 
equation 

\Hcos9)\ = ^. (106) 



where h(cos9) is defined in eq. (p9|). 

It is interesting to note that, in absence of LOFF pairing, i.e. if the 
pairing region vanishes and cos# a = cos#b, one gets for the Meissner and 
Debye mass the same results obtained in the hard dense loop approximation 
for the model with two flavors in the normal (not superconducting) state, 
see e.g. ||27|| . On the other hand if the blocking region vanishes (cos# a = — 1, 
cos 9b = +1) one is in the 2SC phase and eq. (|105|) gives mo = as expected. 



Next consider SH 1 ^ (p). The contribution of the blocking region vanishes 
as discussed in Appendix 3. Therefore the only non vanishing contribution 
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to SH^{p) comes from the pairing region. It could be computed exactly, 
but the corresponding expressions are cumbersome and shed little light on 
the underlying physics. Therefore we prefer to work in the approximation of 
small momenta (\p\ <C A). We find 

-5W ab {p) = 5 ab —^ ^ - 2 



That is 



v; pairing ejj 



2 2 

' v; pairing eff 



(107) 



g 2 [i 2 R 2 f d cos 8 d<fi vtv 

' 3A% 4 J pairing 8^ (<Jfl[/l(C0S 0)]) 



PiPj (108) 



where v = (sin 9 cos <f), sin 9 sin <ft, cos 9) and n = q/q along the z— axis. The 
integration domain is defined by cos9 a < cos 9 < cos Of,. We get therefore 

-*n2(p) = W (/(i2)pi + ^)pf) , (109) 

where k is given by (|94| ) and 



/(i?) = - / dcos#(l -cos 2 #) . (110) 



pairing 



g (R) = - f dcos9 cos 2 9 (111) 

2 J pairing 

are functions of the parameter i? and are reported in fig. 3 (q ~ 1.2|5/i|). 

It is interesting to note the anisotropy of the polarization tensor exhibited 
by these results. One has always g > /; for large R, and neglecting Sfx/fj, 
corrections, one finds approximately 

g(R) 2 



q_ 

M, 



(112) 



Let us finally write down the remaining components of the polarization 
tensor. From ( 107 ) we get 



.2, ,2 



n£(p) = sj£ E VJ ^vl = kp 2 (f(R)(S il 6 jl +5 i2 S j2 )+g(R)5 i3 5 j3 

v; pairing 

(113) 
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Fig. 3 Plots of the functions f(R) and g(R). 

and 

- n«(p) = k PoI P (f(R)(S il 5 jl + 5 i2 6 j2 )+g(R)5 i3 5 j3 ) . (114) 

These results complete the analysis of the LOFF model in the one plane wave 
approximation. From we get the dispersion law for the gluons at small 
momenta. The lagrangian at one loop isQ 

£ = —F^F^-^AlAl (115) 

(sum over the repeated color indices a,b = 1,2,3). Introducing the fields 
Ef = F£ and Bf = ie ijk F" k , and using (jT09|) , ( [TTBD and ([114]) we can rewrite 
the lagrangian ( |115| ) as follows 

£ = \ (ey ^ - 5? 5?) + \rnlA° a Al , (116) 

where 



fl + kf(R) 

l + kf(R) j . (117) 

V l + kg{R) 

This means that the medium has a non-isotropic dielectric tensor e and a 
magnetic permeability A = 1. These results have been obtained taking the 

2 We do not include here the 3 and 4-gluon vertices that however can be handled as in 
[p6[, with the result that the local gauge invariance of the one- loop lagrangian is satisfied. 



25 



total momentum of the Cooper pairs along the z direction. Therefore we 
distinguish the dielectric constant along the z axis, which is 

e,| = l + kg(R) , (118) 

and the dielectric constant in the plane perpendicular to the z axis 

e ± = l + kf(R). (119) 

This means that the gluon speed in the medium depends on the direction of 
propagation of the gluon; along the z axis the gluon velocity is 

V\\ - , (120) 



while for gluons which propagate in the x — y plane we have 

v± ~ . 1 (121) 
f kf(M 



and in the limit of large R, and neglecting 5fi/fi corrections, 

U|| — > —j= tan 9 q v± . (122) 
v2 

with cos^q defined in eq. fl30l) . The anisotropy in the gluon velocities has its 
pendant in the anisotropy of the phonon velocity pi] . 



6 Gluon dynamics in the LOFF phase: cubic 
structure 

The condensate in this case is given by eq. (pJ), so that we will use the 
results of sec. ^| with A e ff given by fl5T|). The calculations are similar to 
the previous case and, similarly, the SU(2) C gluons have vanishing Meissner 
mass and exhibit partial Debye screening. However the dispersion law of the 
gluons is different. 

As a matter of fact we write the one loop lagrangian for the SU(2) C gluons 

as 

£ = l -{EtEt-B«B«) + 5C, (123) 
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with 



(124) 



In the approximation \p\ <C A, SU. 1 ^ is again given by eq. (|107|) , but now 
A e ff is given by (|51~1). One gets 



5C = E?E b , 5, 



: • - , ~ab „ 9 
07T 2 



d cos 9 ( 



VjVi 



2, ,2 



pairing 



in 



'A 2 



+ A a A h 6, 



ab ' 



l eff 



An 2 . / blocking 



dcos8 



Evaluating the integrals one finds 



C 



EtBf) + -M D 2 A%A% 



with the tensor e 1 - 7 given by 



Sy [1 + kt(R)} 



and 



M, 



cos ^ a — cos 6b 



D 




(125) 



(126) 



(127) 



(128) 



where cos 6^ are solutions of eq. ( |106| ). It should also be noted that the 
values of the parameter R for the cube and the plane wave can be different. 

The isotropy of these results can be easily understood noticing that SC 
is a symmetric function of its arguments. As a matter of fact it must be 
symmetric under the residual symmetry Oh and therefore it must be built 
up from the three elementary symmetric functions ^(Xi), I±{Xj) and Ie(Xi), 
defined in eq. (|66|), with Xi = Ef. The invariance of SC follows from the 
invariance of these functions under the three elementary rotations Ri and the 
inversion /. However it is clear from (|125|) that SC is second-order in E l a and 
therefore the only function that can be involved is 12- In conclusion in the 
cubic lattice one necessarily gets 



e ij ~ Oij — A5ij 



(129) 



From eqs. (|125 ) and ([51]) one has that the function t(R) is given by the 
formula 



t(R)5i 



3R 2 
Arc 3 



dcos 



ViVj 



pairing 



2 • 



(130) 
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t(R) can be computed by observing that, using the properties of the 5r 
function 



3R 2 J-y r , „ , , ViVj 
t(R)bij ps — / dcos9d(p rr - — — ^ i2 

47T , J pairing 



k=l J Paring [Sr^V • fljfe))] 

E + [g(R) - / Wl^nf ) • (131) 



fc=l 

Using (RSI) we finally get 



t(R) = -[2f(R)+g(R)] . (132) 



A plot of the function t(R) is in Fig. 4. 




Fig. 4 Plot of the function t(R). 

We finally note that, even if the crystalline structure is not isotropic, in 
eq. ( |126|) the 0(3) invariance is restored. This happens because the tensor 
iij is proportional to the identity. Therefore the dielectric properties of the 
medium will be isotropic and the velocity of propagation of the gluons will 
be the same in all the directions. We stress that these results are valid in the 
approximation of soft gluon momenta. 
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7 Conclusions 



We have used the High Density Effective Theory (HDET) formalism to cal- 
culate the low energy properties of phonons and gluons in the LOFF phase 
of QCD.The HDET formalism is very useful to simplify the dynamics as it 
allows direct elimination of the degrees of freedom which are far from the 
Fermi surface. It bears useful resemblances with the Heavy Quark Effective 
Theory of general use when dealing with heavy hadrons. We have first made 
the modification for the case of a single plane wave condensate and we have 
derived the quark propagators. After that we have gone to the relevant case 
of the cubic strucure of ref. [TT| corresponding to the sum of eight plane 
waves. We have added the suitable term in the Lagrangian bearing the sum 
over the eight plane waves of the cubic structure and derived for this more 
complex case the quark propagators to be used in the calculations. 

The three phonons correspond to the local deformations from the fluc- 
tuations of the crystal as defined by the cubic condensate. They are the 
Nambu-Goldstone bosons from the breaking of translational symmetry. We 
have studied the low energy limit, where wavelengths are much longer than 
lattice spacing, and built up an effective lagrangian for the phonons, which, 
at the second order in the derivatives, depends on three parameters, in agree- 
ment with the general theory of the elastic waves in crystals [24]. 

The parameters of the effective phonon lagrangian are calculated by in- 
troducing external fields with the phonon quantum numbers and performing 
the derivative expansion of the generating functional. The phonons remain 
massless, as expected. The final results for the parameters are remarkably 
simple, the only quantity appearing in the result being the ratio of 6/jl to the 
common length of the eight wave vectors appearing in the condensate. 

We have also derived the effective lagrangian of the gluons of the unbro- 
ken SU(2) C subgroup. The gluon self energy diagrams at one loop can be 
evaluated from the derived effective action. The Meissner mass vanishes as 
in the case of complete isotropy; on the other hand there is a partial De- 
bye screening originating from the existence of blocking regions on the Fermi 
spheres. For soft gluon momenta, in the case of a single plane wave, the 
dielectric constant of the medium is non isotropic and it takes on different 
values according to its direction; also the gluon velocity is anisotropic. In 
the case of the cubic crystal the gluon dispersion law for small momenta is 
however isotropic. It is interesting to remark that at least for gluon long 
wavelengths, the essential anisotropy of the cubic structure is washed out, 
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while it remains visible for the plane wave crystal structure. 
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Appendix 1: The cube symmetry group Oh 

The cube symmetry group Oh is the finite group containing the transforma- 
tions that leave the cube (see Fig. 1) invariant. It is defined as follows: 

O h = O x a t (133) 

where O is the octahedron group (see below) and Cj = {E, 1} is the center of 
symmetry of the group, containing the neutral element E and the inversion 
J: 

-10 \ 

-1 , (134) 

, - 1 / 
i.e. 

E : X\ — > xi, x 2 — > x 2 , x 3 — > x 3 (135) 
I : x\—> —xi, x 2 —> —x 2 , x 3 — > —x 3 . (136) 

The group Oh contains 48 elements, comprising the 24 elements of the octa- 
hedron group O and other 24 given by g ■ I (g G O). The octahedron group 
is defined by the following system of axes of symmetry: 

• Three 4-fold axes^, denoted as C4 in Fig. 1; 

• four 3-fold axes, denoted as C 3 in Fig. 1; 

• six 2-fold axes, denoted as C 2 in Fig. 1. 
Its 24 elements are divided into five classes: 



3 For any crystal a n-fold axis is such that the crystal assumes the same position after 
a rotation of an angle 27r/n around it. Such a rotation is denoted as C ra ; for extension we 
denote by C n not only the rotations but also the corresponding symmetry axes. 
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• Class I, containing only E; 

• class II, containing three rotations C4 and three rotations Cf around 
the three 4- fold axes 4 ; 

• class III containing three rotations Cf three rotations of 2 x 2tt/4 = n 
around the three 4-fold axes; 

• class IV, containing four rotations C3 and four rotations Cf around the 
four 3-fold axes; 

• class V, containing six rotations C 2 around the six 2-fold axes; 

It is easy to construct all the 24 group elements. Class II comprises the three 
C 4 rotations, i.e. rotations of n/2 around the three 4-fold axis: 

/0 - 

#2=01 i . n, 
\i 

and the three Cf rotations: 

C4 : Rl, i?2, R-l (138) 
Class III comprises the three Cj rotations: 

C4 : Rf, i?2, -R3 (139) 
Class IV comprises four C 3 rotations: 

C 3 : A^R^, A 2 = R 1 A 1 R~\ A 3 = R 2 A 1 R^\ A 4 = A 3 A X A^ 

(140) 

and four Cf rotations: 

C3 : A\ , A\ , A\ , A4; (141) 
finally class V comprises six C 2 rotations: 

C 2 : = AiR 2 , 5 2 = A^xAr 1 , B 3 = A 2 BxA 2 x , 

4 In general C™ denotes a rotation of an angle m x 2jr/n around a n-fold axis; clearly 
C" = 1. 
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S 4 = A 3 B 1 A~ 1 , B 5 = A^Al 1 , B 6 = R 1 B 1 R± 1 . (142) 



We finally observe that from the general theorem for finite groups assert- 
ing that the sum of the squares of the dimension of the representations, rij, 
must be equal to the order of the group (in this case 24): 

E«h 24 (143) 

i 

the only possible solution for O is 

ni=n 2 = 1, n 3 = 2, n 4 = n 5 = 3. (144) 



Appendix 2 

In this Appendix we evaluate the coefficients appearing in eq. 



(fc)^(fc) 
j 



k=l 



dv 
8^ 



vgv m 8 iv ■ n k 



5/j, 



To this end we make use of the result 



ViVj 5 I cos 6 — — 1 = f — ViVj 5 | cos 6 — — 
, 3 \ q I J 8tt 3 I q 



S/j, 
Q 



sm 



cos 2 8„ 



— {SnSji + 5 i2 5 j2 ) H — -S i3 S j3 



(145) 



(146) 



with 



We therefore obtain 



cos 9„ 



tfi _ v + 3cos 2 g g -l \ 

P^m - 2^ "^V Rqm^ij ^ g °pq H g ^p3Cg3 J , 



(147) 



(148) 



where the eight matrices .R^ are rotation matrices that transform the eight 
vectors into the unit vector (0, 0, 1). They are given by 



( 


1 


1 

V2 





\ 




1 


1 


[2 








f 


V 3 




\ 






1 











(149) 
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and, in the notations of the Appendix 1 



= rWr-\ rW = rWrI, r^ = rVr 3 , 

R^ = rMb 3 , R® = R®B4, R<n=R<$B 3 , R® = R^B 4 .(150) 



On the other hand the eight matrices C\ are 



C 



(k) 



( i vi k) 4 k) \ 



(k) 



(k) 



(k) (k) , 

V m Vi 1 J 



(151) 



with 



(k) 

m 
(k) 

V2 

(k) 



;+i,-i,-i,+i,-i,+i,+i,-i) 
;+i,+i,-i,-i,-i,-i,+i,+i) 
;+i,-i, +1,-1, +1,-1, +i,-i) 



Let us evaluate f3 lm . We notice that 



C 



(k) 



<% + E 

r=l,2,3 



where 
and 



P'i 



SisSjt + SuSj S , r, s, t in cyclic order 



E ^ fc) = o, E « fe) = 85 r 

fe=l,-,8 



fe=l,-,; 



(152) 

(153) 

(154) 
(155) 



Using the orthogonality of the R^ matrices we can evaluate the first term 
of the r.h.s. of ( |148| ) immediately 



aii _ c ; n 2 flA r , 3 cos2 B q ~ 1 y> p (fc) p (fc)^(fc) 

° k=l 



(156) 



From the very definition of the matrices i?^ we get 



D (fc) (fc) 



Ji3 



dW ( fc ) j 



(157) 
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Therefore 

= sin 2 e q S lm S t3 + 3CQS2 ^~ 1 E <7«<7« . (158) 

Z4 fc=i 

Using the expression ( |153| ) for C^- we find 

/C = f W« + 3C ° S2 3 ^" 1 E P£nP& • (159) 

Appendix 3 

For the case of a single plane wave we get the following expression for the 
polarization tensor: 



- n£(p) = z5 afe ^ E [ w» + v^/iCp) + (vv" + v v v») Up) 

V 

- iSa^^J^-V^V-VW), (160) 



where 



e ,V-lV-(l + p) 



[ ^ = J ie WT^dW ' (161) 
4 M = / ft JM5M ' (163) 

. - + « 2 /;; (2 , +t ^, )gM . d«) 

and -D(£) = V -£V •£ — A^ + ze. We note that the contribution proportional 
to J in ( |160| ) arises from the diagram of fig 2 b and corresponds to the 
lagrangian £ 2 m ©■ 

The integrals I\ (p) and Ji (p) are infrared divergent in the blocking region 
where A e ^j = 0. One has to regularize them going to finite temperature and 



imaginary frequencies (£q — » 2mT(n + 1/2)) (for a similar case see ||20|| , eqns. 
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(2.193) and (2.194)). The results that we get for the blocking and the pairing 
region are respectively: 

ii(0) = h{p) = Ji(0) = hip) = -2ttz, J 2 (0) = hip) = 0, J = —fi 2 (165) 
and 

h(0) = h(p) = h(0) = h(p) = -m, l 2 (0) = -m, J = y^- ( 166 ) 
From these we can easily compute the Debye and Meissner masses using 

ti; blocking 



+ 2g tiU -V"V U -VY" (167) 



,2, ,2 



- n^(o) . . = £ [vv + t>v" 

pairing 271 - — L 

pairing 

+ 2g^ - 2 (V»V U + V U V^) ] . (168) 

We get therefore the vanishing of the Meissner mass 

ng(0) = 0. (169) 

On the other hand we get a non-vanishing Debye screening mass from the 
blocking region: 

- Il2(0) = 2W«/£ £ =5 ab ^4[ dcosO , (170) 

7T -. , . Z7T Jblockinq 

v; blocking M 

where 

r rcos9 a rl 

/ dcos6= / dcos6> + / rfcos6». (171) 

J blocking J — 1 J cos Of, 

The domain of integration of the blocking region has been obtained by the 
condition 

R\h{cos6)\ > | , (172) 

which gives, using the approximate expression for h(cos8) (0), and for R > 
tt/2: 

cos9 q . J cos 9 q . _ 

COStf a = COS 6^ = 171171 < 1, ^- f , (.173) 

1+ 25 I 1 ~ 2i? 
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where cos^g = Sfi/q. From the equations ( |170| ) and ( |172| ) one gets the Debye 
mass: 

gfi f cos0 a -cos0 6 
m v = — + 2 • (7) 

These results can be easily generalized to the case of the cube. 
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